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& ' Abstract 

h- » : 

qv ■ A new formula for the nuclear charge radius is proposed, dependent on the mass 

number A and neutron excess N — Z in the nucleus. It is simple and it reproduces all 
the experimentally available mean square radii and their isotopic shifts of even-even 
nuclei much better than other frequently used relations. 
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1 Introduction 

The aim of the present investigation is to find a phenomenological macroscopic formula 
for the nuclear radius. A charge distribution inside a nucleus depends not only on it's 
mass number A and the ground state deformation caused by shell effects, but also on the 
neutron excess N — Z. The equilibrium deformation of a nucleus can be established by 
the minimalization of the potential energy obtained by one of the microscopic methods. 
But to know the range (R) of nuclear charge in any (i?, ip) direction, it is necessary to 
describe properly the radius (i?oo) of the corresponding spherical nucleus which has been 
deformed keeping its volume constant. 

The simplest known liquid drop formula for the nuclear radius 

Ro = r ■ A 1 / 3 (1) 

with tq = 1.2 fm doesn't lead to the experimentally known mean square radii (r 2 ) (MSR) 
[1] or their isotopic shifts [2]. A popular liquid droplet model [3] gives also rather unsat- 
isfactory results for MSR. 

We are going to find the best dependence of R 00 on mass number (A) and neutron 
excess (N — Z). A corresponding formula which reproduces the experimentally measured 
isotopic shifts of MRS was already proposed in refs. [4, 5]. That formula was obtained 
without the fit to the absolute values of the model independent of MSR data [1]. Now 
we would like to propose a new simple expression for the proton radius (i?oo) which 
reproduces not only the isotopic shifts but also the magnitude of MSR for all even-even 
nuclei with proton numbers Z > 8. 

2 Description of the calculation 

The MSR values of even-even nuclei are calculated macroscopically with a uniform 
charge distribution 

J |s fo r r < R f . 

P °-{ for r>R (2) 

where V is the volume of a nucleus conserved during deformation 

4 „ r 2w r l r R({i3 x y,v,<p) 

V=-nBL= dip dcostf r 2 dr . (3) 

3 Jo J-i Jo 

R is the radius dependent on deformation, calculated with the equilibrium shapes {/3a} 
taken from ref. [6]. 

The MSR values are given by the integral 

/■2ff ,1 rR({p x }-A<p) 

(r ) = po / dp d cos •& / r dr (4) 

Jo J-i Jo 

and they are evaluated for axially symmetric nuclei with the radius 



R({P X },$)=R ({[3 X }) 
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l + £/3 A P A (costf) 

A=l 



(5) 



where Pa are the Legender polynomials and 



MiPx}) = Roo L J^ 



l + 2> A P A (cos#) 



A=l 



-1/3 



dcosfiy . (6) 



The last relation originates from the volume conservation condition when a nucleus de- 
forms and i?oo is the radius of the corresponding spherical nucleus. For practical use the 
formula (4) for small quadrupole deformations could be estimated approximatly by: 



(r2) = ^of 1+ 8 y + 
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where Qoo — | e -^-Roo * s the electric monopole moment of the spherical nucleus and Q20 
is the electric quadrupole moment of the deformed nucleus. 

The charge radius i?oo m formula (6) is chosen in the following form (see also [4-5]) 

k N - Z^ 



Roo = r A^ (l + - - a— — J 



The set of r , K and a parameters are obtained by the fit of (r 2 ) ex P [1] and its isotopic 
shifts 5 < r 2 >^ [2] for even-even nuclei with the proton numbers Z > 8. We have 
assumed that the sum of the mean square errors: 
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should be minimal. Simultaneously we have kept the minimal sum of square errors of 
M SR isotopic shifts for all experimentally measured sets of even-even isotopes 

5(r 2 ) A ' A ' = (r 2 ) A - (rY (10) 

i.e. 
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The results obtained with our simple formula for the charge radius (8) will be compared 
in the next section with the liquid droplet estimates of MSR [3] . 



3 Results 

After an analysis of 22 isotopic shifts of MSR for even-even nuclei we have proposed in 
Ref. [5] a formula for i?oo containing the isotopic term, 

R 00 = r A 1 ^ (l - «^-^) (12) 

with the parameters a = 0.2 and r = 1.25 fm, but without the ~ 1/A term, i.e. k, = 
0. This formula was good enough to reproduce all the MSR isotopic shifts of even- 
even elements above Z > 38 and the average magnitude of MSR for the gold isotopes. 
Unfortunately it was rather impossible to get a reasonable agreement of MSR values of 
lighter elements 8 < Z < 38 without introducing an additional k/A term. Fitting R 00 
described by formula (6) to the MSR [1] and 5{r 2 ) A ' A [2] data for all the Z > 8 elements 
we have recorded the following set of parameters: 

r = 1.240/m, K = 1.646, a = 0.191 

. In every case the sum of square errors of MSR is more than six times lower than in the 
liquid droplet model [3]. It is interesting, that the a parameter is almost stable around 
a ~ 0.2, independently on the number of fitted data and inclusion of the ~ 1/A term. 
The values r and a are almost the same for the case of Z > 38 elements without the k/A 
in R 00 . 

We have also tested the different powers of A in the new term k/A u (n = |, |, • • • 2). 
We have achives the best agreement for n — 1, but for n = | and | the quality of the 
fitment was comparable with that for n — 1. 

The results of our investigation are illustrated in the figures 1-5. In the Figs. 1-3 
the results of the MSR values obtained in the three cases of i?oo are illustrated. 

In Fig. 1 the MSR values obtained with the traditional liquid drop formula (1) 
are plotted. They are obviously to small even for the spherical case (dashed line) in 
comparison with the experimental data [1] marked by stars. The difference between the 
macroscopic results including the effect of the deformations of nuclei taken from Ref. [6] 
and the experimental data (r 2 ) exp are represented by the length of vertical lines shown 
under each experimentally known element ("bars"). The differences are rather large, 
especially for the light elements. 

In Fig. 2 the MSR results were obtained using our previous formula [5], but with 
the better set of parameters (r = 1.256/m, a = 0.202), fitting, both < r 2 > and 
5{r 2 ) A,A values. The fitment is much better than that in Fig. 1, both for heavier and 
lighter elements, even in the spherical case. When including deformation, the differences 
(bars on the figure) between macroscopic and experimental value (r 2 ) t h — {r 2 ) e xp almost 
disappear, but the lighter elements are still not well reproduced. It was necessary to 
introduce an additional term in the formula for the charge radius in order to increase the 
model estimates of MSR for the lightest nuclei. 

In Fig. 3 the MSR results for the formula (8) with the best set of parameters fitted 
to the (r 2 ) exp and 5{r 2 ) A ^ A data are shown. The agreement is good. In spite of a few 



cases the differences between the macroscopic and experimental MSR values have almost 
disappeared when the deformations of the nuclides were taken into account [6] . 

One can compare these results with the liquid droplet calculation according to ref. [3] 
illustrated in Fig. 4, where the theoretical results are much smaller than the experimental 
data. Also the results for MSR obtained with the newest version of the droplet model [6] 
are only a little better than those of [3] in Fig. 4. It is obvious that a new analysis of 
liquid droplet parameters in order to reproduce MSR values is needed. We hope that our 
investigation will be helpful. 

The importance of introducing (TV — Z) dependence in i?oo is illustrated in Figs. 
5-6, where the MSR isotopic shifts S(r 2 ) A ^ A of all the experimentally known even-even 
isotopes [2] are compared with our estimates. 

In Fig. 5 the differences between the macroscopic results calculated with equilibrium 
deformations of [7] and experimental data [1] 

A(r 2 ) A > A '=5(r 2 )t h A '-S(r 2 ) A £ 

are plotted in the case of the traditional liquid drop formula (1). The differences are 
large. They reach even 1 fm 2 for some nuclei. When including the i?oo depending on 
isospin ^-^- (Fig. 6), these differences decrease significantly. One has to add that the 
same quality of the fit of 5{r 2 ) A ' A as with formula (8) was already achieved for the two 
parametric formula (12), but the magnitude of MSR values for light nuclei are not well 
reproduced. 

4 Conclusions 

We can draw the following conclusions from our investigations: 

• The traditional liquid drop formula, depending on mass number A only, is not able 
to give a good mean square nuclear charge radius (r 2 ) or its isotopic shifts. 

• In order to reproduce the data of (r 2 ) and 5{r 2 ) of the even-even nuclei one should 
use the three parametric formulas for the spherical nuclear charge radius 

„,,, / 1.646 N-Z\ „ 

i?oo = 1-240 A 1 '* (l + —j- + 0.191 — -^ — J fm 

• For the nuclei with Z > 38 one can omit the ~ 1/A term and use the formula 

Roo = 1.256 A 1/3 (l - 0.202 — -— ) fm 

• The liquid droplet model gives systematically larger MSR estimates of the nuclear 
charge. 
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6 Figures Captions 

Fig. 1. The mean square charge radii in fm 2 of even even nuclei obtained without 
deformation (dashed line) with the traditional formula liquid drop compared with ex- 
perimental data [1] (crosses). The differences between the macroscopic results evaluated 
including equilibrium deformations [6] and the experimental data [1] are illustrated by 
the bars around the ordinates line. The macroscopic MSR values are to small. 

Fig. 2. The same as in Fig. 1 but with the isospin dependence in the i?oo formula 
(12). The differences for light nuclei are still of the order ~ lfm 2 . 

Fig. 3. The same as in Fig. 1 but with our new formula for R 00 (8). The 
experimental data is almost ideally reproduced. 

Fig. 4. The same as in Fig. 1 for the liquid droplet macroscopic results with 
parameters of Ref. [3]. The theoretical values are to small up to 2fm 2 . 

Fig. 5. The differences between the liquid drop and the experimental isotopic shifts 
of the MSR of the charge in fm 2 . The discrepancy between these macroscopic results 
and the experimental data exceed lfm 2 . 

Fig. 6. The same as in Fig. 5 for the three parametric i?oo formula (7). The 
deviations between our theoretical estimates and experimental data do not exceed 0.5 fm 2 . 
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